function sMatrix = construct_sparse_graph(graphNodes,graphEdges)
tStart = tic;

p = [];
for i = 1 : length(graphEdges(:,1)) * 2
    p = [p; 1];
end

sMatrix = sparse([graphEdges(:,2); graphEdges(:,1)], [graphEdges(:,1); graphEdges(:,2)], ...
    p, length(graphNodes(:,1)), length(graphNodes(:,1)));

tElapsedCMatrixBuilding = toc(tStart)